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Abstract 

The drag and diffusion coefficients of charm and bottom quarks propagating through quark 
gluon plasma have been evaluated for conditions relevant to nuclear collisions at Large Hadron 
Collider and Relativistic Heavy Ion Collider. Both the radiative and coUisional processes of 
energy loss are included in evaluating the effective drag and diffusion coefRcients. The Landau- 
Pomeronchuk-Migdal (LPM) and the dead cone effects on the radiative energy loss of heavy 
quarks have been included. With the effective transport coefficients the Fokker Plank (FP) 
equation has been solved for the evolution of heavy quarks in QGP. The solution of the FP 
equation has been used to evaluate the nuclear suppression factor, Raa and the elliptic flow, 
for the non-photonic single electron spectra resulting from the semi-leptonic decays of 
hadrons containing charm and bottom quarks. It is observed that the experimental data from 
RHIC on Raa and can be reproduced simultaneously within the pQCD framework for the 
same set of inputs. 
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1. Introduction 

The nuclear collisions at Relativistic Heavy Ion Collider (RHIC) and the Large Hadron 
Collider(LHC) energies are aimed at creating a phase where the bulk properties of the 
matter are governed by the (hght) quarks and gluons. Such a phase of matter is called 
Quark Gluon Plasma (QGP). The study of the bulk properties of QGP is a field of great 
contemporary interest and the heavy flavors, mainly, charm and bottom quarks, play a 
crucial role in such studies. As the relaxation time for heavy quarks are larger than the 
corresponding quantities for light partons, the light quarks and the gluons get thermalized 
faster than the heavy quarks. Therefore, the propagation of heavy quarks through QGP 
may be treated as the interactions between equilibrium and non-equilibrium degrees of 
freedom and the Fokker-Planck (FP) equation provides an appropriate framework [1] 
for such studies. In this work we would like to evaluate the elliptic flow and nuclear 
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suppression factor of heavy flavours within the framework of FP equation. The results 
of the calculations will be compared with the available experimental data. 

The evolution of heavy quarks momentum distribution function, while propagating 
through the QGP are assumed to be governed by the FP equation, which reads. 



where the kernels Ai and Bij are given by 



(1) 



Ai = j d^kuj{p, k)ki 

Bij= j (fikco{p,k)kikj. (2) 

iov \ p \^ 0, Ai ^ jPi and Bij DSij, where 7 and D stand for drag and diffusion 
co-efficients respectively. 

The basic inputs required for solving the FP equation are the dissipation co-efficients 
and initial momentum distributions of the heavy quarks. The drag and diffusion coef- 
ficients have been evaluated by taking in to account both the coUisional and radiative 
processes [2]. In the radiative process the dead cone and LPM effects are included [2]. 
In evaluating the drag co-efficient we have used temperature dependent strong coupling, 
as{T) [3]. The Debye mass, ^ 9{T)T also a temperature dependent quantity used as a 
cut-off to shield the infrared divergences arising due to the exchange of massless gluons. 
The initial momentum distribution of heavy quarks in pp collisions have been taken from 
the NLO MNR [4] code. The solution of the FP equation for the heavy quarks is con- 
voluted with the fragmentation functions of the heavy quarks to obtain the transverse 
momentum {pr) distribution of the D and B mesons. For heavy-quark fragmentation 
function, the Peterson function has been used. The pT distribution of the electrons from 
the semi-leptonic decays of D and B meson are evaluated using the standard techniques 
available in the literature. The ratio of the pt distribution of the electron from the decays 
of heavy flavours produced in heavy ion collisions to the corresponding (appropriately 
scaled by the number of collisions) quantities from the pp collisions is defined as the 
nuclear suppression factor, R^a'- 

coll d^pj-dy 

which will be unity in the absence of re-scattering. The STAR [5] and the PHENIX [6] 
collaborations have measured the Raa for non-photonic single electron as a function of 
Pt for Au-I-Au at y/s^N = 200 GeV. 

The experimental data from both the collaborations [5,6] show Raa < 1 for pt > 2 
GeV indicating substantial interaction of the heavy quarks with the plasma particles. The 
spectra evaluated using the formalism described above reproduces the data reasonably 
well (Fig. 1). It is clear from the results displayed in Fig. 1 that the charm quarks suffer 
more suppression than bottom quarks at RHIC energy. 

Now we discuss the elliptic flow resulting from non-central collisions of nuclei. When a 
heavy quark propagates along the major axis of an ellipsoidal domain of QGP (resulting 
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Fig. 1. The variation of Raa sjid withpy. The theoretical results are obtained within the frame- 
work of FP equation described in the text with the same set of inputs. The drag and diffusion coefficients 

have been calculated using pQCD approach. Contributions from charm and bottom quarks arc shown 
separately. The initial temperature and thermalization time are taken as 400 MeV and 0.2 fm/c respec- 
tively. 



from the non-central collisions) then the number of interactions it encounters or in other 
words the amount of energy it dissipates or the amount of momentum degradation that 
takes place is different from when it propagates along the minor axis. Therefore, the mo- 
mentum distribution of electrons originating from the decays of heavy flavoured hadrons 
produced from the fragmentation of heavy quarks propagating through an anisotropic 
domain will reflect such anisotropy. Thus the spatial anisotropy due to non-central col- 
lisions will be reflected through the momentum distributions. The degree of momentum 
anisotropy will depend on both the spatial anisotropy and more importantly on the cou- 
pling strength of the interactions between the heavy quarks and the QGP. The drag and 
diffusion coefficients depend on the temperature of the background medium (QGP) which 
evolves in space and time due to expansion. Therefore, the drag and diffusion will also 
change due to the flow of the background. The flow of the background has been treated 
within the ambit of (2+1) dimensional hydrodynamics [7]. The coefficient of elliptic flow, 
is defined as: 

(Pt) = {cos{2(j))) = r J dN I (4) 

J "'f' dydpTd4,\v=0 

We evaluate in the current formalism [8] (see also [9]) and compare the results with 
experimental data [10] (Fig. 1, right panel). 

The prediction for the nuclear suppression factor and elliptic flow of the heavy quarks 
to be measured at LHC energy through the semi-leptonic decays are depicted in Fig. 2. 
The sensitivity of the results on the equation of state (EoS), i.e. on the velocity of sound 
is also considered. Lowering of Cg gives more suppressions as observed in Fig. 2. Lower 
value of velocity of sound, makes the expansion of the plasma slower enabling the 
propagating heavy quarks to spend more time to interact in the medium and hence lose 
more energy before exiting from the plasma which results in less particle production at 
high Pt- 
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Fig. 2. The variation of Raa and v^^ with px for LHC for 0-10% centrahty. The initial temperature 
and thermalization time are taken as 700 MeV and 0.08 fm/c respectively. The theoretical results are 
obtained within the framework of FP equation described in the text with the same set of inputs. The 
drag and diffusion coefficients have been calculated using pQCD approach. 

Some comments on the -Raa vis-a-vis are in order here. The Baa contains the 
ratio of pt distribution of electron resulting from Au-|-Au to p+p collisions, where the 
mimerator contains the interaction of the heavy quarks with the medium (QGP) and 
such interactions are absence in the denominator. Whereas for both the numerator 
and the denominator contain the interactions with the medium, resulting in some sort of 
cancellation. Therefore, the K factor dependence for the two cases seems to b{^ different. 

In summary, we have evaluated the nuclear suppression, Raa and the aziniuthal asym- 
metries {v2^) using the Fokker-Planck equation. The transport coefficients are evaluated 
using pQCD interaction. We have shown that our results on both the Raa and can 
reproduce the experimental data with the same set of inputs. 
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